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Hot-hole lasers in III-V semiconductors. 
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Department of Applied Physics, Technical University Delft, Lorentzweg 1, 2628 CJ DELFT, The Netherlands. 

Following the success of p-Ge hot-hole lasers, there is also potential for using other semiconductor materials, 
notably III-V's such as GaAs and InSb. Previous analysis had suggested that a large effective mass ratio between 
the heavy and light holes is advantageous, which implies that InSb would make an excellent hot-hole laser. Using 
our Monte Carlo simulations of both GaAs and InSb hot-hole lasers in combination with a rate equation model, 
we see that previously accepted criteria used to predict performance are not always reliable, and we suggest 
suitable alternatives. The simulation results include gain and gain bandwidth as a function of field strength and 
laser frequency, and alternative field orientations and photon polarizations are considered. Comparisons are 
made with bulk p-Ge systems. The optimum conditions predicted by our simulations could then be used in the 
design of quantum-well hot-hole lasers. 



This paper was published as J. Appl. 
(2001). 



I. INTRODUCTION 



Phys. 90, 1692 



Hot-hole lasers emit in the THz (far-infrared) with an un- 
usually broad gain spectrum, allowing amplification and gen- 
eration of laser pulses on a picosecond time scale [Jj], 0, 51 
B> § @> tyP e °^ l aser nas been realised in bulk 

p-doped Germanium (p-Ge), and produced gains of about 
0.25cm _1 . Since the THz band has important potential ap- 
plications in medical imaging and short-range (e.g. office) 
communications, and due to the widespread applications of 
GaAs H [Io[ [llj, there has been interest in the potential of III- 
V materials such as GaAs and InSb as hot-hole lasers. Most 
predictions have been based on simple scattering rate and ef- 
fective mass arguments (see e.g. [|I]); so to investigate the 
possibilities of III-V materials more thoroughly we have done 
a set of Monte Carlo simulations for both GaAs and InSb. We 
cover a range of field strengths, orientations, and doping con- 
centrations. 

Hot-hole lasers usually consist of a p-doped Ge crystal in a 
cryostat cooled to about 20K, with crossed electric and mag- 
netic fields applied, and are described extensively in Ref. |Jy]. 
In this system, the ideal lasing cycle for a hole is depicted in 
Fig. |l]as follows: (1) the electric field accelerates a heavy hole 
to above the optical phonon energy Elo, (2) the hole scatters 
into the light hole band by emitting an optical phonon, (3) due 
to its lighter mass, the hole is now localised on a closed cy- 
clotron orbit in the light hole band, and this localisation forms 
an inversion, and so (4) stimulated emission of a photon trans- 
fers the hole from the light back to the heavy hole band. In 
order to get the streaming motion (1) combined with the cy- 
clotron orbits of (3) it is necessary to apply the correct ratio 
of electric to magnetic field strengths as determined by the 
heavy to light hole effective mass ratio. If this is not done we 



may get cyclotron orbits which do not reach Elo in the heavy 
hole band, or light holes which reach Elo and scatter instead 
of being localised in cyclotron orbits. 

The results we present in this paper are obtained using a 
Monte Carlo simulation (see section ||) which can deal with 
for bandstructure and scattering processes properly. However, 
in order to describe the basic mechanism we now present a 
simple rate equation model. 

The holes are, roughly speaking, in one of three different 
locations: "H" the lower part of the heavy hole band; "H+" 
the upper part of the heavy hole band (above Eio)', and "L" on 
cyclotron orbits in the light hole band. Holes in "H" stream up 
to "H+" with rate r stream jl2[], or their streaming is disrupted 
by intraband scattering with rate rnn- Holes in "H+" will 
most likely emit an optical phonon, and either scatter into "L" 
with rate yrop, or back into "H" with rate (1 — y) rrjp. Holes in 
"L" can either scatter by phonon or impurity processes back 
into "H", or emit a photon and scatter into "H" also. This 
is diagrammatically represented on fig. [l](b), and the popu- 
lations (Nh, Nh+, Nl) in these locations can be described by 
rate equations: 



d 

-rN H = 
dt 

d 

dt N " + = 

d AT 

dt NL = 



-r stream N H + [(1 - i)r 0P + r HH }N H + + r LH Ni\) 
J rr stream N H - [rop + run] N H + , (2) 
+jropN H + - r LH N L , (3) 
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where y is the fraction of optical phonon emissions which 
end up in the LH band. Using a density of states ratio argu- 
ment and ntftH 3> mw> this is {muH / 'tnui)~~ • Note that 
this model says nothing about the inversion, which is localised 
in k space. 

None of the steps in the ideal lasing cycle (1-4) are very 
efficient, and so it is best to regard the hole transport as being 
predominantly a sea of unhelpful scatterings with the lasing 
cycle superimposed on top of it. It is only if we can make the 
lasing cycle strong enough by selecting a suitable material, 
and adjusting the field strengths, field directions, or doping 
concentration will we see lasing. Past analysis of the lasing 
cycle has provided number of criteria that are usually used 
to determine whether steps (1-4) will be sufficiently strong 
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FIG. 1: The hot-hole laser: (a) shows the streaming motion up the 
heavy hole band (1), the scattering to the light hole band (2), the cy- 
clotron orbits in the light hole band (3), and the lasing transition (4); 
(b) shows a simple three-state model of the operation of the laser, 
with the transition rates being: r stream the rate of streaming up HH 
to above the optical phonon energy, yqp the optical phonon emis- 
sion rate, run the HH intra-band scattering rate, and ruj the rate of 
depopulation of the LH band. 



JlH , [l4], 15 1 . In particular, criteria for the temperature and 
electric and magnetic field strengths were given to ensure that 
the heavy and light holes would be unlikely to scatter below 
the optical phonon energy; that heavy holes do scatter when 
above this energy; and that some of these scattered holes reach 
the light hole band. They also favoured a high effective mass 
ratio ntHH/mrji, such as that in InSb, making it easier to en- 
sure that the heavy holes are accelerated to the optical phonon 
energy whilst the light holes remain in closed orbits. We show 
here that such arguments, especially the latter, do not always 
hold true. 



II. SIMULATIONS 

The simulations were done in two stages. First, for 
chosen electric (E) and magnetic (B) field directions (i.e. 



E[0lI]B[100]; £[112]fl[110]; £[2lT]fi[011]; £[lll]fi[112] - 
see fllq l ), one Monte Carlo simulation was done for each 
point of a grid of field strengths covering the range of interest. 
The simulations produced distribution functions of heavy and 
light holes and the average scattering times for the system un- 
der the chosen conditions. The second stage used the distribu- 
tions to calculate inversion-gain cross-sections <j/; and the av- 
erage scattering times were used in a Drude model to calculate 
the optical absorption a a . These two were then combined and 
multiplied by the impurity concentration n, to give the net gain 
per unit length X = (a/ — a )n, as a function of photon fre- 
quency for each grid point. This three dimensional dataset of 
net gains was then analyzed to produce contour plots of maxi- 
mum gain, enabling us to see the best directions and strengths 
of E and B to choose to get optimum gain. 

Our Monte Carlo simulations JT^] o f this system include a 
full k.p band structure calculation ||18|] and all important scat- 
tering processes: optical phonon scattering due to deforma- 
tion potential (OPD) and polar interactions (OPP), acoustic 
phonon scattering due to deformation potential (ACD) and 
piezoelectric interactions (ACP), and (Coulomb) ionized im- 
purity (IIM) scattering ^ BJ. We treat the effect of the elec- 
tric and magnetic fields classically, giving the holes continu- 
ous trajectories in fc-space - comparison with experiment for 
p-Ge systems indicates that this approximate treatment is ad- 
equate for the field strengths we consider. All simulations are 
done at a lattice temperature of 20K. They follow the progress 
of a single hole through a large number of scatterings (typi- 
cally ~ 16000), with the ergodic theorem being used to justify 
the use of the simulation's time-average as an ensemble aver- 
age. Since both the hole-hole scattering processes and ionized 
impurity scatterings are mediated by screened coulomb po- 
tentials, this means that hole-hole scattering can be allowed 
for by doubling the impurity concentration in the simulation - 
since the hole concentration is equal to the impurity concen- 
tration [Q. 

The Monte Carlo simulations used the standard overesti- 
mation technique where for each scattering process the post- 
scattering direction of the hole was chosen at random, and the 
differential scattering rate was overestimated by an isotropic 
rate just higher than its maximum value. However, this proved 
to be very slow for the ionised impurity scattering, especially 
at low impurity concentrations. Therefore, in this case we 
weighted the choice of scattering angle by the angular de- 
pendence of the scattering rate, and thus avoided generating 
a large proportion of inefficient overestimations. 

Optical absorption was estimated using a Drude free car- 
rier model, which uses zone centre effective masses and the 
hole scattering times calculated in the Monte Carlo simula- 
tions. This model assumes that each scattering completely 
randomises the phase of the particle. However, IIM scatter- 
ing, which makes a significant contribution to the average 
hole scattering times, involves a large proportion of scatter- 
ings which only slightly deflect the hole. As a result, the or- 
dinary treatment of IIM scattering would significantly over- 
estimate the scattering rate used to calculate the optical ab- 
sorption. We avoid this problem because our IIM algorithm 
weights the choice of angle, while the scattering rate itself 
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is isotropic. By choosing an algorithm compatible with our 
model of optical absorption, we achieved both an efficient 
model of IIM scattering whilst avoiding the necessity of a 
more complex absorption model. 



III. GERMANIUM 

To validate our model and simulation code, we first did 
some simulations of hot-hole lasers in p-Ge and compared our 
results against those typically given in the literature. The opti- 
mum gain contour plots are shown on fig. |[ and table | shows 
the relative frequencies of the different scattering types at the 
point of optimum gain. The columns ACD, IIM, and OPD on 
this table refer to the acoustic phonon, ionised impurity, and 
deformation potential optical phonon scattering processes re- 
spectively. These simulations also provided us with a good 
reference with which our GaAs and InSb simulation results 
can be to compared. The results, for an impurity concentra- 
tion of Hi = 0.025 x 10 16 cm~ 3 , are shown in Fig. 0, where 
the optimum gain is about 0.25/cm at F 
B = 4T [Oil]. 
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TABLE I: The relative frequencies of different scattering types in Ge 
atapoint of optimum gain: £ = 5kV/cm [2lI],B = 4T [Oil] and«,= 
0.025 x 10 16 . There were a total of 16000 scatterings, and the hole 
was in the LH band for 250 ± AAps, and the HH band for 1076 ± 49 
ps. The average time between collisions in the LH and HH bands 
were tin = 4.9ps and tnH = Lips respectively. NB: mfjH/ m LH '■ 

^stream ~ 1 .3ps 



-1. 



Transition ACD IIM OPD Total 

HH^HH 749 161 13619 14529 

HH^LH 48 8 612 676 

LH^LH 4 95 21 120 

LH^HH 158 90 428 676 



These simulations also enabled us to test the accuracy of 
our simple model of optical absorption against the more so- 
phisticated model of Lok [ |Tr5[ [22], which included second 
order transitions of phonons and photons, or ionised impu- 
rities and photons; while all scattering processes where taken 
into account. We obtained optical absorptions within ~ 25% 
of Lok at photon energies of 5meV - due to the inverse de- 
pendence of absorption with frequency, the fit was much bet- 
ter above 5meV, but worse below. Lok's model has optical 
phonon scattering dominating optical absorption, with IIM 
scattering appearing to have a negligible effect. In his results, 
the lower gains at higher impurity concentrations are due to 
IIM scattering from the light hole to the heavy hole band, 
causing a loss of inversion. This has a more significant effect 
on the optical gain than the extra optical absorption caused 
by the IIM scattering. Using our algorithm that weights the 
choice of IIM scattering angle, we see the same. 
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FIG. 2: Gain in p-Ge, n; = 0.025 x 10 cm , laser photon polar- 
ization in the E direction, (a) The net gain after free-carrier optical 
absorption has been taken into account. The thick contour is at zero, 
and contours step up in gain by 0.025cm~' . (b) The net gain as a 
function of magnetic field and photon frequency (for the optimum 
electric field). 
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IV. GALLIUM ARSENIDE 

The obvious III-V material to investigate is GaAs, as it is 
already widely used for many other applications |lO[ It 
has zone center effective masses of roughly mhh ~ 0.50 and 
mui w 0.07, with a ratio rather similar to that in germanium 
(Ge). The main difference compared to Ge is that GaAs has 
additional scattering processes - it is a polar material, so there 
are also piezo-electric phonon (ACP) processes and polar op- 
tical phonon (OPP) processes to consider. While the piezo- 
electric processes are usually not very significant, the polar 
optical phonons have a significant role, largely because their 
strength is inversely dependent on the phonon momentum. 



TABLE II: The relative frequencies of different scattering types in 
GaAs at a point of optimum gain: E = 8kV/cm [211], B = 7T [01 1] 
and in = 0.025 x 10 16 . There were a total of 16000 scatterings, 
and the hole was in the LH band for 180 ± 25 ps, and the HH band 
for 687 ± 5ps. The average time between collisions in the LH and 
HH bands were t^u = 3.8ps and = 0.7ps respectively. NB: 



m HH I m LH ~ 8; r str . 



1.7ps" 



Transition ACD ACP IIM OPD OPP Total 



HH^HH 899 

HH -> LH 43 

LH -> LH 3 

LH^HH 222 



186 383 4198 9001 14667 

4 34 188 297 566 

7 181 1 9 201 

41 193 8 102 566 



Simulations were done for a range of electric and magnetic 
field magnitudes, from 2kV/cm to 20kV/cm, and 2T to 12T. 
The maximum fields were chosen somewhat higher than is 
usually practicable in experiment to ensure we covered the 
full range of interest. For «, = 0.025 x 10 16 cm~ 3 , the E/B 
ratio for optimum gain was 1 . 14 (cf 1 .25 for Ge). 

A full array of simulations was done for n, = 0.025 x 
10 16 cm~ 3 , an electric field direction of J? [211], and a mag- 
netic field direction of fi[011] (relative to the crystal axes). 
This moderate impurity concentration showed a best net gain 
of 0.08cm _1 at high fields (^8kV/cm, 7T) (as shown in Fig. 
U). Although reasonable inversion-gain occurs down to fields 
of a few kV/cm (or T), the optical absorption is strong enough 
to spoil it except at higher fields. The net gain for a range 
of different field orientations were also checked at the opti- 
mum 8kV/cm, 7T point, but £ [211] &B[011] were confirmed 
as giving the greatest gain. In general, the orientation depen- 
dence can be significant, with sometimes up to a factor of two 
difference in the net gain. 

Next, simulations were done at higher and lower impurity 
concentrations for the best-gain value of E and B fields. The 
higher impurity concentration (n; = 0. 100 x 10 16 cm~ 3 ) results 
indicated that there is at best a marginal gain. This is because 
the increased impurity scattering has two effects. First, the ex- 
tra scattering increases the calculated optical absorption; and 
secondly, impurity scattering is more efficient at depopulat- 
ing the LH band by producing LH — > HH scattering than it is 
at repopulating it with HH — > LH scatterings. A table of the 
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FIG. 3: Gain in p-GaAs, rif = 0.025 x 10 16 cm~ 3 , laser photon polar- 
ization in the E direction, (a) The net gain after free-carrier optical 
absorption has been taken into account. The thick contour is at zero, 
and contours step up in gain by 0.025cm _1 . (b) The net gain as a 
function of magnetic field and photon frequency (for the optimum 
electric field). 
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relative frequencies of the different scatteringtypes at E = 8 
kV/cm [211], B = 7 T [Oil] is given in Table jj| 

The lower impurity concentration simulations (n, = 0.005 x 
10 16 cm~ 3 ) show a better net-gain cross-section 07 — 0,4. How- 
ever, since the number of holes in the simulation is related 
to the impurity (doping) concentration, fewer impurities also 
means fewer holes, and fewer holes means that the gain per 
length X is proportionally reduced. So although the net-gain 
cross-section is up to double that for n; = 0.025 x 10 16 cm~ 3 , 
the net-gain per cm reduces to X ~ 0.04cm _1 . 

The dependence of the net gain on the polarization of the 
output light is small. Typically one orientation (of B, E, or 
E x B) is better than the other two, which are otherwise very 
similar. Output light polarized in the B direction requires the 
hot-hole laser to be in the Voigt configuration, but E polarized 
output can be achieved by either the Voigt or Faraday configu- 
rations. In our simulations, the polarization in the electric field 
(E) direction gave the best gain; this contrasts with the com- 
ments in Ref. [|l3[ saying that a polarization parallel to the 
magnetic field (B) minimises heavy hole scattering and hence 
the optical absorption. 



V. INDIUM ANTIMONIDE 

An often mentioned candidate for a III-V hot-hole laser is 
p-InSb, because its effective mass ratio mHu/mLH is much 
greater than that in either GaAs or Ge. The simple argument 
is that since the light hole cyclotron orbits can be more tightly 
confined in fc-space, this leads to a more pronounced inver- 
sion. However, since our simulation code cannot treat light 
hole Landau levels, we restricted our simulations to magnetic 
fields under 2T; and so were unable to make predictions at the 
higher fields where the confined light hole cyclotron orbits 
might give a significantly increased gain. 



TABLE III: The relative frequencies of different scattering types in 
InSb at a point of optimum gain: E= 1.8kV/cm [011], B= IT [100] 
and m = 0.025 x 10 16 . There were a total of 16000 scatterings, 
and the hole was in the LH band for 55 ± lips, and the HH band 
for 1326 ± 13 ps. The average time between collisions in the LH 
and HH bands were tm = 4.5ps and tnH = 1.4ps respectively. NB: 
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Transition ACD ACP IIM OPD OPP Total 
HH -> HH 1681 47 648 3610 9679 15665 
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From Fig. [| we can see that the peak gain at n, = 0.025 x 
10 16 cm~ 3 is about 0.05cm _1 , and the trend is the usual one 
for better inversion at higher fields. The dependence of the 
net gain on the polarization of the output light is small, as for 
GaAs, and indeed Ge. For InSb, the best field orientation is 
not an electric field direction of [211], and a magnetic field 
direction of [011]; but £[011] B [100]. A table of the relative 
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FIG. 4: Gain in p-InSb, n, = 0.025 x 10 16 cm 3 , laser photon polar- 
ization in the E direction, (a) The net gain after free-carrier optical 
absorption has been taken into account. The thick contour is at zero, 
and contours step up in gain by 0.025cm _1 . (b) The net gain as a 
function of magnetic field and photon frequency (for the optimum 
electric field). 
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frequencies of the different scattering types at E = 1 .8 kV/cm 
[Oil], B = 1.0 T [100] is given in Table g. 

As with GaAs, simulations were done for InSb at the same 
higher and lower impurity concentrations («,■ = 0. 100, 0.005 x 
10 16 cm~ 3 ) for the best-gain value of electric and magnetic 
fields. The higher impurity concentration produced a predic- 
tion of no net gain - for the same reasons as when the higher 
impurity concentration reduced the gain in GaAs. The low 
impurity concentration gave reduced gain, with the resulting 
smaller hole concentration again overcoming the slight in- 
crease in gain cross section. 



VI. DISCUSSION 
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as important (GaAs), or was regarded as bad (InSb). This is 
in contrast to simple arguments which state that large ratio is 
good because it helps ensure efficient streaming up the HH 
band whilst retaining good cyclotron orbits in the LH band - 
unfortunately, too large a ratio brings the LH density of states 
down too low, inhibiting population transfer to the LH band 
and hence inhibiting inversion. The simple arguments also 
miss the influential role of the dielectric constant in control- 
ling impurity scattering rates (GaAs). 

Next, we need to consider the criteria given in past analyses 
as described in section |. Here we can clarify their roles using 
the simple rate-equation model (Eqs. [j]-^j) described earlier. 
For Nl <C Nh +Nh+, the fraction of population in the LH 
band in the steady state is 



These results indicate a number of important considerations 
for hot-hole laser operation in bulk III-V materials. Firstly, 
control of the impurity concentration is critical, as e.g. in 
GaAs, the material can go from no net-gain through its peak 
back to no net-gain with only a factor of 20 or so range in 
impurity concentration. Secondly, as in Ge, the orientation of 
the electric and magnetic fields with respect to the crystallo- 
graphic axes is important - while the difference in the amount 
of net gain from the best choice in InSb to the second best 
choice is not large, the region over which it can be obtained 
is. This then is particularly important if there are restrictions 
on the maximum fields allowed. Finally, it is clear that the 
performance in p-GaAs and p-InSb is not as good as in p-Ge 
hot-hole lasers; the main reason being that the non-polar p- 
Ge does not have the polar optical phonon (OPP) and acoustic 
pizeoelectric phonon (ACP) scattering processes, and that this 
extra scattering has a more significant negative effect than the 
benefits of the improved effective mass ratio. 

For GaAs, OPP scattering transfers more holes into the 
LH band than deformation optical phonon (OPD) scattering 
- however, it also scatters more out. As a result, the net 
transfer of holes into the LH band is similar for both opti- 
cal phonon processes. The much greater LH — > HH rate for 
OPP is due to the small momentum transfer Ak in such pro- 
cesses: the OPP rate is inversely proportional to Ak. Also 
note that the GaAs ionised impurity scattering (IIM) rates are 
roughly double those for Ge - this is largely due to the differ- 
ence in dielectric constants, which appear to the fourth power 
in the denominator of the scattering rate prefactor. The di- 
electric constant for Ge is about 16, compared to 13 for GaAs. 
Also, the dielectric constant affects the screening length which 
also alters the IIM scattering rates. Unfortuately it is not pos- 
sible to easily compare the IIM scattering rates and dielec- 
tric constants between Ge or GaAs and InSb, as their effec- 
tive mass ratios are so different - (niHH / Q e q u ^ s ~ 8, but 
(m H H/m LH ) !nSb w 30. 

In InSb the net transfer of holes into the LH band is sim- 
ilar for both optical phonon processes. Also the HH — > LH 
scattering rates in InSb are consistently lower than for GaAs, 
which is due largely to the lower LH density of states caused 
by InSb's lower LH effective mass. 

Note that for both GaAs and InSb in the above discussion, a 
large effective mass ratio yiihh I m LH was either not mentioned 



N H + Nh+ r LH (r p + r stream + r HH ) ' 

The three important criteria are: (i) a large y, which means 
a small rn.Hn I m LH ratio (but not so small as to disrupt either 
the HH streaming or LH cyclotron orbits); (ii) that thh should 
be as small as possible compared to r stream + rop', and (iii) that 
should be much less than the smaller of r stream or rop. We 
have extracted these scattering rates from our optimum gain 
Monte Carlo simulations, and using these in eqn.(Q) gives the 
same trends as is observed for the net gain - P = 0.10, 0.08, 
0.05 for the Ge, GaAs, and InSb simulations respectively. 
Also, the simulation y's were generally close to the ratio of 
effective masses. 

In summary, we have presented a rate equation model of 
hot-hole lasers and simulation results that clarify the roles of 
performance criteria presented in earlier analyses. We show 
that in general III-V materials should have an acceptable per- 
formance, which contrasts with previous criteria that sug- 
gested more optimistically that (e.g.) InSb should make an 
excellent hot-hole laser. This demonstrates the need to treat 
the lasing cycle as a whole, as in our rate equation model, 
rather than treating each step in isolation. Finally, the perfor- 
mance of III-V hot-hole lasers might most easily be improved 
by moving to a modulation doped quantum wells - thus elim- 
inating the impurity scattering by moving the impurities away 
from the active region, and giving us opportunity to engineer 
the band structure. 
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